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1. Introduction 



The Skyrme model is a three-dimensional field theory with topological soliton so- 
lutions, which is believed to provide a good model of nuclear physics; the topolog- 
ical solitons (skyrmions) are identified with baryons. In this paper we investigate 
skyrmion chains, in other words topological solutions which are periodic (or quasi- 
periodic) in one of the spatial directions. 

Isolated skyrmions have been well-studied, up to relatively high charge ^. 
Multiply-periodic configurations have also been considered, for example a "Skyrme 
domain wall" solution (doubly-periodic) f^, and a "Skyrme crystal" (triply-periodic) 
[0, H, 1^. However, to date no one has considered topologically-nontrivial chains. 
There have been studies of cylindrically-symmetric skyrme strings P, 0, but these 
are topologically trivial. 

The field equations of the Skyrme model are non-linear, and no exact analytic 
solutions are known; so the only direct means of studying solutions is by numerical 
simulation. Since the system is three-dimensional, numerical methods are not fast. 
However, a number of alternative methods are known: these typically involve a 
relatively simple approximation, which is empirically observed to approximate the 
true (numerically-determined) behaviour of the solutions. The three most prominent 
ansatze of this type are the product ansatz (see for example 0), the Atiyah-Manton 
construction [^, and the rational map ansatz [|^. We have been able to adapt the 
first two to study chains, but the rational map ansatz appears unsuited to this task. 

The plan of the paper is as follows. Section 2 contains a discussion of the topology 
of skyrmion chains; section 3 has results on adapting the Atiyah-Manton construction 
to obtain approximate skyrmion chains from Yang- Mills calorons; section 4 discusses 
the approximation of skyrmion chains in terms of parallel vortex-antivortex pairs; 
and section 5 contains some full 3-dimensional numerical results. 

We show that there is a 1-skyrmion chain with a preferred period (the period 
which minimizes the energy-per-period of the chain), and that the caloron-generated 
field is a good approximation to this solution. The vortex-antivortex pair is not 
as good, but becomes more accurate for denser (lower-period) chains. Generally 
speaking, a chain of 1-skyrmions is unstable for larger periods, in the sense that the 
individual 1-skyrmions tend to clump into higher-charge skyrmions. We give some 
numerical results which illustrate this, in particular showing the appearance of chains 
of 2-skyrmions and of 4-skyrmions. 
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2. The Skyrme model and chains 

The static Skyrme field U is defined on M.^, and takes values in SU(2). Defining 
Li = U^^dU/dx\ we take the energy density of U to be 

£ := -^Tr(L,L,) - lTr([L„ L,][L„ L,]); (2.1) 

and the normalized energy of U on M.^ is 

E := / Sdx^dx^dx^. (2.2) 

The boundary condition at spatial infinity is f/ ^ 1 as r ^ oo, where 1 denotes 
the identity element of SU(2). So topologically, a Skyrme field defines a map S^ -^ 
SU{2), and such a map has a degree B E Z. This topological charge can be computed 
by the integral 

B= Bdx^dx^dx^, (2.3) 



where 

S = ^e,,fcTr(L,L,Lfc) (2.4) 



is the topological charge density. The energy ( p.2| ) satisfies the topological lower 
bound E > B. Finite-energy fields U which locally minimize E will be called 
skyrmions. The symmetry group of the system (including the boundary condition) 
consists of the spatial translations, the spatial rotations 0(3)sp, and the iso-rotations 

0(3)iso. 

A Skyrme chain with period j3 > and relative orientation R G S0(3)iso is a 
finite-energy Skyrme field U : M.^ ^ SU{2) satisfying 

U{x,y,z + (3)=R-U{x,y,z), (2.5) 

minimizing the energy functional 

E = / £dx dy dz, 

Jo Jr2 

and satisfying the boundary condition U —>■ 1 as x"^ + y'^ ^>- oo. One may think of 
such a configuration as a chain of equally-spaced skyrmions along the ^-axis, with 
each skyrmion being iso-rotated by R relative to its neighbours. In the Appendix, it 
is shown that the quantity 



B = / Bdxdydz 

'o Jr2 
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is necessarily an integer, and this integer is called the topological charge of the chain. 
(If R is the identity, then it is already clear that B is an integer, since then the field 
U is in effect a map from S"^ x S^ to SU(2), and B is the degree of this map.) 

A single skyrmion has the spherically-symmetric hedgehog form f/(x) = exp {if{r)x^a^}, 
where a^ denotes the Pauli matrices. Consequently, the iso-rotations R G S0(3)iso of 
a single skyrmion may be identified with the spatial rotations, and we may think of 
the relative orientation of two skyrmions as being a relative orientation in physical 
space. If we have two well-separated 1-skyrmions, then it is well-known p| that the 
force between them depends on this relative orientation: 

• If the two skyrmions have the same orientation, then they repel. This case is 
called aligned. 

• The strongest repulsive force occurs when one of the skyrmions is rotated by 
vr about the line / joining them. This case is called the repulsive channel. 

• The strongest attractive force occurs when one of the skyrmions is rotated by 
vr about an axis perpendicular to /. This case is called the attractive channel. 

The symmetry of such a 2-skyrmion configuration can be inferred from the dipole 
model of skyrmions [0, or equivalently by looking at a superposition of hedgehog 



configurations. A pair of skyrmions which are aligned, or in the repulsive channel, 
has an axial symmetry 0(2) (consisting of rotations and reflections which fix /); but 
a pair in the attractive channel has only a discrete symmetry group D2 (generated 
by reflections in two perpendicular planes whose intersection is I). 

Guided by this, we define three types of unit-charge skyrmion chain as follows. 
Writing f/(x) = exp{w(x)}, where w takes values in the Lie algebra su{2), we may 
regard the relative orientation operator R defined in ( |2.5| ) as acting on w(x) via the 
adjoint action. Suppose that along the z-axis, w has the form w = g{z)v, where v is 
a fixed element of the Lie algebra su{2), and g{z) is a real- valued function. We then 
identify the following three special cases: 

• A chain with i? = 1 is called aligned. 

• A chain for which i? is a rotation by n about v G su{2), is called maximally- 
repulsive. 

• A chain for which i? is a rotation by vr about an axis perpendicular to f G su{2) 
is called maximally- attractive. 
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As before, the maximally-attractive chain has only a discrete symmetry, whereas the 
other two have a continuous axial symmetry. We expect that only the maximally- 
attractive chain will have an energy less than that of an isolated skyrmion, and it is 
this type that we shall concentrate on in what follows. 

3. Skyrme chains from calorons 

The Atiyah-Manton ansatz is method of generating approximate skyrmion configu- 
rations on M.^ by evaluating the holonomy of 4-dimensional Yang-Mills instantons. 
The construction is topologically natural, in the sense that the holonomy of an A^- 
instanton is an A^-skyrmion; and for an appropriate choice of the instanton scale, the 
resulting skyrmion field is a surprisingly good approximation to the actual skyrmion. 
The method involves evaluating the holonomy of the gauge field along a family 
of parallel lines in the 4-dimensional space. For example, we could choose the family 
of lines parallel to the x°-axis in M^, and then the Skyrme field is obtained via the 
path-ordered exponential integral 



\J —oo 



tAj ] (J/iA/ 



In one periodic version of this construction which has been studied previously, it 
was shown that one may obtain a good approximation to the Skyrme crystal from 
instantons on the 4-torus ||T2| . 

If we want to obtain a Skyrme chain, then we should start with a gauge field on 
S^ X M? satisfying 

i\n,\X ~r P •; X , X , X ) ^ ±L ' J\n[X , X , X , X ), ( O.l 1 

and integrate along a family of parallel lines perpendicular to the x^'-axis. Calorons 
provide examples of gauge fields satisfying (|3.1|). Let A^ be the gauge potential of a 
caloron, strictly-periodic with period /?, with instanton charge equal to 1, and with 
vanishing monopole charge (see |13|, Q for details). So near infinity, Aq does not 
wind, and has the form Aq th ifia^, where < fi < n/ 13. Now make the gauge 
transformation A'^ = gA^g^^ — {dug)g^^, where g = exp(iyux''cr^). In the new gauge, 
Aq — !► at infinity, and the gauge field is no longer strictly-periodic, but satisfies 



A'^(x° + p, x\ x^ x^) = hA^{x°, x\ x^ x^)h 



5 



where h = exp(iyLi/3cr^). This gauge choice is known as the algebraic gauge. We 
can obtain 1-skyrmion chain configurations by computing the holonomies of such 



-5- 



p 


12n^E2 


1271^^4 


B-1 


0.2 


10.56 


509.2 


2.0 X 10-^ 


0.3 


16.61 


314.3 


1.6 X 10"^ 


0.4 


23.20 


213.7 


1.1 X 10-^ 


0.5 


29.68 


160.4 


7.7 X 10-5 


0.6 


35.40 


135.4 


6.3 X 10-5 


0.7 


40.21 


124.6 


9.3 X 10-5 


0.8 


44.25 


119.6 


1.9 X 10-^ 



Table 1: Caloron approximation 



calorons in the algebraic gauge. Analytic expressions are known for all calorons 
with unit instanton charge and vanishing monopole charge |T3]; so it is feasible to 
study Skyrme chains using this family. Since we are mainly interested in maximally- 
attractive chains, we restrict attention to the calorons which yield R^ = 1 but R y^ 1; 
these are the calorons with fi = 7r/2/3. By contrast, calorons with yU = 0, such as the 
Harrington- Shepard calorons [|I6|, give rise to aligned chains. 

The calorons we are interested in are symmetric under rotations about an axis 
in M^. If we choose to evaluate holonomies along lines parallel to this axis, the 
Skyrme chains will also have an SO {2) symmetry: these will be maximally-repulsive 
chains. If, on the other hand, we evaluate holonomies along lines perpendicular to 
the symmetry axis, we obtain maximally-attractive chains. 

We have implemented the Atiyah-Manton construction for maximally-attractive 
chains numerically, and evaluated the energies of the resulting skyrmion chains. The 
family of calorons we used is parametrized by a scale parameter p and a period (3. 
The caloron with scale p and period (3 is in fact a rescaling of the caloron with scale 
p/(3 and period 1, and since the components of the Skyrme energy behave simply 
under rescalings, it was sufficient only to consider calorons with fixed /? = 1 and a 
range of values of p. 

The holonomies were evaluated using the Runge-Kutta method. We evaluated 
energies in a finite box —L < x,y < L, and extrapolated in both the box size and 
the lattice spacing to obtain energies accurate to within 0.1%. We also calculated B 
to check the accuracy of our method. Our results are summarized in Table 1. 

In figure 1(a), we have plotted the minimum energy of this approximate skyrmion 
chain, as a function of the period (3. The graph was obtained by interpolating the 
data in Table 1 to obtain E2 and E4 as polynomial functions of p, and then minimizing 
the energy E = E2J3 + -E4//? with respect to variation in p. In particular, we see that 
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the energy of these caloron- derived configurations has its lowest value ii^ ^ 1.16 for 
period (3 ^ 2.1. 
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Figure 1: The energy of a 1-skyrmion chain, versus the period (3, for the caloron con- 
struction (solid curve), the vortex ansatz (dashed curve), and from a full 3D calculation 
(dot-dash curve). 



4. The vortex ansatz 



In this section we describe an alternative ansatz for chains, which is based on the 
idea that a chain, especially for small period, splits into constituents. We define a 
Skyrme vortex to be a field of the form 



Uy = exp [^{9 - iyz)ia^] exp [/(r)icr^] exp [^{6 + i^z)icr'^] , 



(4.1) 



where {r,6) are polar coordinates in M^, and z/ is a positive constant. The profile 
function /(r) is required to satisfy the boundary conditions /(O) = 7r/2 and /(r) -^ 
as r ^ oo. Note that U^ is smooth on M.^, and periodic in z with period 2-11 /v. The 
energy density of U^ is 



£=('f')2 , ^/ , ,,2.;„2^ , /^/n2COsV , / x/N 2. .2 „:„ 2 x , ^^^ slu V COsV 



+ z.^inV + (/'r 



+ (/')VsinV + 



■ (4.2) 



The profile function /(r) is chosen to solve the Euler- Lagrange equation for j^ Sr dr. 
But note, for example from the second term in S, that U^ has infinite energy per 
unit period. This is a consequence of the fact that U^ has the non-constant form 



Uy ^ exp{i6a^) 



(4.3) 
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as r — >■ CX3. The Skyrme vortex (4J.) has been used before, in a different context [17 
constructing vortex loops by taking finite lengths of Skyrme vortex and joining their 
ends together. 

Notice that the field U{x,y, z) = f/t,(x, —y, —z) winds in the opposite direction 
to Uy at infinity. So we can obtain a field which is constant at infinity by taking a 
superposition of two vortices via the product ansatz 

U = U^U2 

with Ui(x, y, z) = f/„(x — a, y, z) and f/2(x, y, z) = Uy{x + a, —y, —z), where a is a 
positive constant. We can obtain a similar, but more symmetric, field by using the 



relativized product ansatz [|18 



U = (f/it/2 + t/2f/i)/v/(det(t/it/2 + f/2f/i). 

The superposition satisfies the boundary conditions of a maximally-attractive chain 
with period P = n/u and R = diag(— 1, —1, 1) G S0(3), and has the same symme- 
tries. The field resembles a pair of parallel vortices separated by a distance 2a, and 
its topological charge is 1. 

When the separation 2a of the vortices is large, they attract each other, as 
the following heuristic argument shows. Let C > be sufficiently large for the 
approximation (|4.3| ) to be valid for r > C, and let a be larger than C. We can 
evaluate the energy of the superposition by splitting M? into three regions: the two 
discs of radius C centred on the vortex locations, and the exterior. The energy 
within each disc tends to a constant as a tends to infinity. The energy in the exterior 
diverges as a ^ oo; a calculation shows that the leading contribution at large a is 
47r/3 In a. So one can reduce the energy by reducing the separation 2a, as claimed. 

So far we have not justified our choice of superposition procedure: it is important 
to ask whether there is another way to superpose two vortices to obtain a lower 
energy. Again, we have a heuristic argument why our superposition is the right 
thing to do, at least for large separation. Let F denote the exterior of the two discs 
Di, D2 of radius C and centres {x,y) = (±a, 0), as before. Let ip : F ^ U{1) be a 
map such that ip\dDi has winding number 1 and ip\dD2 has winding number —1. We 
want to minimize the energy 

e = p [ ||^-^#|| Vx. 
Jf 

The ansatz used above corresponds to taking ip = exp[i(6'i — 6*2)], where 9i{x, y) is the 
angle between (x — a, y) and the x-axis, and 6^2 (x, y) is the angle between [x + a, y) 
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Table 2: Vortex approximation 



and the x-axis. If the energy of this field is close to the true minimum, then we know 
our ansatz is a good one. Notice that the Skyrme term has disappeared from our 
energy functional; this is because the Skyrme term evaluates to zero for any f/(l) 
field. The easiest way to find the minimum energy is to stereographically project 
from M^ to S*^; the energy e is conformally-invariant, so we are allowed to do this. 
The stereographic projection can be chosen so that the two circles are described 
by ^ = a and ^ = vr — a in spherical coordinates 9 G [0,7r], G [0,27r), where 
sin(Q;) = C /a. The energy functional is now written as 



/9 



27r 



{{r'dei'f 



sm 



-2 



e{i>-^d^i)f)sm9d9d(t). 



A Bogomolny argument shows that this energy is minimized by ilj{9,(f)) = exp(i0); 
and the minimum energy is 47r/91ncot(Q;/2). For large a, this agrees with our super- 
position, to leading order. 

We have evaluated the energy of the superposition of two vortices for a range of 
values of j3 and a. The energies were evaluated in a finite box, and we extrapolated 
in the box size and the lattice spacing to obtain results accurate to within 0.1%. We 
tried using both the product ansatz and the relativized product ansatz, and found 
that the energies obtained agreed. We also evaluated the topological charge B as 
a check on our methods. Table 2 shows the minimum energy of the superposition, 
together with the value of a for which this energy is attained. The energy is plotted, 
as a function of the period /3, in Figure 1(a). 
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5. Numerical Minimization. 

In order to assess the approximation schemes of the previous two sections, and also 
to investigate the stabihty of the maximally-attracting 1-skyrmion chain, we imple- 
mented a full 3-dimensional numerical minimization of the energy E. The method 
involved a first-order finite-difference scheme for E, with the spatial points {x, y, z) 
being represented by a rectangular lattice having lattice spacing h. The boundary 
condition U = 1 was imposed at |x| = L, \y\ = L. The energy was minimized using 
a conjugate-gradient method. The errors in ii^ as a result of the finite lattice spacing 
h and finite size L go like h"^ and 1/L^ respectively, and we obtained the /i — > 0, 
L ^ oo limits by extrapolation. The resulting values for E have an error less than 
0.05%. 

The first simulation looked at the maximally-attractive 1-skyrmion chain, for 
each of the periods f3 = 1.9,2.0,2.05,2.1 (since all the indications are that the pre- 
ferred period /3min is close to 2). The results are plotted in Figure 1(b), together 
with a parabola fitted to the resulting four points; we see that the preferred period 
is /3min = 1.98, and that the minimal energy is -Emin = 1.143. 

For large period, we would expect the maximally-attractive 1-skyrmion chain to 
be unstable to clumping, for the same reason that a finite collection of separated 
1-skyrmions will clump together. For small period, however, the 1-skyrmion chain 
might be stable. To investigate this, we first looked at the periodic 2-skyrmion 
chain, where the initial configuration was taken to be (a deformation of) a pair 
of skyrmions on the 2;-axis, in the attractive channel. As before, the numerical 
simulation involved fiowing down the energy gradient, to reach a local minimum. 
The results are illustrated in Figure 2, for two values of the period, namely [3 = 4 
and f3 = 5. The 1-skyrmion chain with period [3 = 2 (essentially the preferred period) 
corresponds to a 2-skyrmion chain with period /3 = 4, and this 2-skyrmion chain turns 
out to be stable. The minimal configuration is depicted in the upper row of Figure 2. 
Subfigure (a) plots the function cr = ^ tr(f/), or rather the surface (t(x) = 0; subfigure 
(b) plots the charge density B, or rather the surface i3(x) = 0.2 x maxB. 

For the larger period [3 = 5, however, the two skyrmions coalesce to form a single 
2-skyrmion, which (as one would expect from the M^ case) has a toroidal shape. This 
is depicted in the lower row of Figure 2, where the same quantities are plotted as 
in the upper row. The instability which leads to this coalescence is rather weak: in 
fact, the energy of the toroidal chain with period [3 = 5 is only slightly (0.1%) less 
than twice the energy of a 1-skyrmion chain with period 2.5. 

Although the 1-skyrmion chain with [3 = 2 is stable when considered over two 
periods, one might expect that it will in turn be unstable to clumping if examined 



-10- 



(a) (jforp=4 



(b) B for p=4 




(c) oforp=5 




• 




M 



• 




Figure 2: The field component a and the charge density B for the 2-skyrniion chains with 
periods /3 = 4 and /? = 5. 



over more than two periods. For example, there is evidence suggesting that a 4- 
skyrmion chain (chain of a-particles) might be particularly favourable: finite-length 
chains of a-particles occur in R^, at least when one modifies the system by adding a 
significant pion mass [^, |2^, ^. We again used numerical minimization, this time 



looking for the minimal-energy chain of 4-skyrmions, with various periods. We use 
the periodicity condition which amounts to rotating neighbouring a-particles by vr 

with i?2 



about the periodic axis |21|; this is of the form (|2.5| ) with R^ = 1 and R j^ 1. The 
a-particle chain with period (3 = 8 has has significantly (4%) lower energy than four 
times the energy of a 1-skyrmion with period 2; the plot of its charge density in 
Figure 3(a) shows that, for this value of the period, the a-particles are quite well 
localized. One can lower its energy by reducing the period, and so allowing the a- 
particles to move closer together: in Figure 3(b), we see the a-particle chain with 
period /3 = 3, and its energy is 6% lower than the minimal value {E = 1.143) for a 
1-skyrmion chain. 



6. Concluding Remarks 

From Figure 1 we see that the Atiyah-Manton construction gives a reasonably good 
approximation to the minimal-energy 1-skyrmion chain, with energy only 1% above 
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(a) B for p=8 

e 



(6) B for p=3 




Figure 3: The charge density B for 4-skyrmion chains with periods /? = 8 and /? = 3, over 
three periods. 



its true value. The vortex-antivortex approximation is not quite as good, but does 
emphasize the constituent structure of skyrmion chains. It appears to be a common 
feature of topological soliton systems that soliton chains exhibit a constituent struc- 
ture when the soliton size is comparable to the period: the solitons fragment into 
fractional-charge objects, and this fragmentation occurs in a transverse direction, so 
that rotational symmetry about the chain axis is lost. These features are apparent 
from Figure 2(b), where one sees evidence of the parallel vortex-antivortex pair. For 
low-period chains (more precisely, where the period is small compared to the natural 
soliton size), one expects the chain to resemble a vortex pair of this type. 

The Skyrme crystal p, ^ Q has an energy-per-baryon oi E = 1.036, and one 
could construct a skyrmion chain by cutting a chain out of this crystal, in other 
words by truncating it in the x- and y-directions. A very thick chain obtained in this 
way should have rather low energy (becoming lower as the chain became thicker). 
Among thin chains (with transverse size comparable to the natural soliton size), 
the chain of a-particles may well be the lowest-energy solution. The preliminary 
calculation reported in the previous section showed that its energy is E' ~ 1.07 for 
period /? = 3. It would be worth making a more comprehensive study of A^-skyrmion 
chains for various A^, for a wide range of periods (3, and with various periodicity 
conditions; but that will require more intensive computational effort than we have 
used in deriving the results reported here. 
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7. Appendix: Topology of Chains 

In this Appendix, we show that the topological charge B is an integer, even when 
the field is not strictly-periodic. It is a special case of the following generalisation of 
the degree theorem, which appears to be new. 

Theorem. Let S be an n-dimensional compact manifold without boundary, with 
volume form uj, and with iJ"~^(S) = and iJ"(S) = Z. Let M be an (n — 1)- 
dimensional manifold such that H^{M x S^) = Z. Suppose that 5*0(2) acts on S 
and that uj is S'0(2)-invariant. Fix an element a G 5*0(2) and let : M x R ^ S 
be a map satisfying 

(f){x, y + 13)= a(t){x, y)\/x e M,y e R. 

Then (p has an integer degree, computed by the integral 

Furthermore, deg{(f)) is independent of the choice of S'0(2)-invariant volume form u. 

Proof. The idea of the proof is simple: we deform to a strictly-periodic map using 
the 5*0(2) action, and show that the integral is unchanged by this deformation. 
First, we introduce some notation: we write the 5*0(2) action as i?s : S ^ S, with 
s G M/Z ^ S0{2). Let X G TS be the associated vector field. Let to e M/Z be 
such that RtQ = a^^, and let t{x,y) = toy/P be a function on M x M. We define a 
deformation 

0(x,|/) = Rt(x,y)i(l)ix,y)), X e M,y eR. 

Then is a strictly periodic map, hence has a degree computed by 

Now we show that 

/ [ ^*cu= [ [ 0*0.. 

Jo JM Jo Jm 

For any form 9 G A*S, one can show that 

0*^ = (j)*Rie + (l)*{ixR*t0) A dt. 
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Here ix denotes the inner derivative of a form. In the particular case 9 = u, one has 
R^uj = uj (because the volume form is S'0(2)-invariant). Hence 

0*u; = 0*u; + (f)*{ix^) A dt. 

By Cartan's formula, we have 

LxU! = ixduj + dixoo, 

where Lx denotes the Lie derivative. We see immediately that duj = 0, because 
UJ G A"S. On the other hand, Lxuj must vanish since uj is S'0(2)-invariant. It 
follows that ixuj is closed. Since iJ"^^(S) = 0, ixuj is exact, in other words, there 
exists a yU G A""^S such that ixi-o = dfj.. Therefore 

0*u; = 0*a; + 0*(rf/i) A dt 
= (f)*uj + d{(f)*fiAdt). 



Integrating and applying Stoke's theorem, we obtain 



/ (t)*UJ = I 0' UJ 

Jm Jo Jm 



(f)*fi A dt 



M 



The boundary term vanishes since t is constant on the domain of integration, so we 
have the desired result. 

That deg{(f)) is independent of the choice of volume form follows from the corre- 
sponding property of the classical degree. If u' is any other S'0(2)-invariant volume 
form, then 

/ (f)*uj' = / / 0V 
Jm Jo Jm 

[ 0*^ 
JM 

/ (l)*UJ. 

Jm 
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